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GEOMETRICALLY INDUCED TWO-PARTICLE 
BINDING IN A WAVE GUIDE 

HELMUT LINDE 


Abstract. For mathematical models of quantum wave guides 
we show that in some situations two interacting particles can be 
trapped more easily than a single particle. In particular, we give 
an example of a wave guide that can not bind a single particle, 
but does have a geometrically induced bound state for two bosons 
that attract each other via a harmonic potential. We also show 
that Neumann boundary conditions are ‘stickier’ for two interact¬ 
ing bosons than for a single one. 


1. Introduction 

Over the last two decades a considerable amount of research has 
been done on mathematical models for quantum wave guides (see e.g. 
man□ei and references therein). Typically a particle in such 
a structure is modelled by a Schrodinger operator on some tube-like 
domain in two or three dimensions. The main object of interest is the 
spectrum of these operators, and especially their low-lying eigenval¬ 
ues which indicate the presence of bound states for the particle. Such 
trapped modes have been proven to exist, e.g., for tubes with local de¬ 
formations, bends, or mixed boundary conditions. Much less is known 
though about the binding of several interacting particles in such set¬ 
tings Tin, HI Edi • In [[I^ Exner and Vugalter addressed the question 
how many fermions can be bound in a curved wave guide if they are 
non-interacting or if they interact via a repulsive electrostatic poten¬ 
tial. It is clear that for these systems a smaller number of particles 
can be bound more easily than a higher number of particles. In the 
present article we consider the somewhat opposite case and show that 
under certain conditions two bosons with an attractive interaction can 
be bound more easily than one particle alone. 

Our work is inspired by the analogous effect for Schrodinger opera¬ 
tors * 1 in free space: Consider for a particle of mass m the operator 

H = -—'—A + V(x) 

2m 

in L 2 (M n ) with a non-trivial, compactly supported and bounded poten¬ 
tial V < 0. It is well known that for n > 2 the attractive potential V 

This work was supported by CONICYT. 

1 We choose units in which the Planck constant K is equal to one. 
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may be too weak to have bound states, i.e., H may not have negative 
eigenvalues. If this is the case, the same potential may still give rise 
to bound states of a system of two particles that attract each other. 
This can be understood by physical intuition if one assumes that the 
two particles act in some sense like one particle of the double mass. 
After all, as far as the existence of eigenvalues is concerned, doubling 
the mass has the same effect as doubling the strength of the potential. 
In the present article we discuss whether an analogous effect can occur 
for purely geometrically induced bound states in wave guides. 

More precisely, we describe a quantum mechanical particle in a wave 
guide by the Dirichlet Laplacian —A in L 2 (h2), where Q is a straight 
strip or tube. The spectrum of this operator is purely continuous and 
contains every real number above some threshold, which is the lowest 
eigenvalue of the Laplace operator on the cross section of Id. It is 
known that geometrical perturbations like bending the tube or local 
deformations of the boundary can give rise to eigenvalues of —A below 
this threshold. In analogy to the case of the Schrodinger operator with 
a weak attractive potential, we ask the following question: Does a wave 
guide exist that doesn’t have a bound state for one particle, but that 
does have a bound state for a system of two interacting particles? 

This question is not so easy to answer by physical intuition, because 
the existence or non-existence of geometrically induced bound states 
for one particle doesn’t depend on the mass of the particle in question. 
This means that the intuitive ‘double mass argument’ for two parti¬ 
cles in an attractive potential doesn’t apply to this situation. Despite 
that, we will show in the following two sections that the answer to the 
question above is ‘yes’ by giving an appropriate example. 

2. Two-particle bound states in deformed wave guides 

We assume our wave guide to be the domain Slcl 2 given by 

^ = {( a ?/): \y\ < / 0 )} 

where 

( | for |x| > L/ 2, 
f(x) — < e/2 fora; = ±L/2, 
y h/ 2 for |x| < L/2 

with L > 0, h > 1 and 0 < e < 1. We impose Dirichlet conditions on 
dtt, which includes the ‘barriers’ at x = ±L/2. Our geometry can be 
interpreted as a cavity of length L and width h coupled weakly (if e 
is small) to two semi-infinite straight wave guides. We choose to set 
rn — |, such that the one-particle Hamiltonian is simply Hi = —A. 
Then standard arguments imply that 

Oess(#l) = [7T 2 ,00). 
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Figure 1 . Sketch of the wave guide 

Eigenvalues may occur depending on the choice of the parameters L, h 
and e, but we will show: 

Lemma 2.1. If L~ 2 + h~ 2 > 1 then for small enough e there are no 
eigenvalues of H 1 below if 1 , i.e., in this case the wave guide has no 
one-particle bound states. 

On the other hand, we consider a system of two bosons of mass 
m — which interact via the harmonic potential 

V = a(x! - x 2 ) 2 + cx{y 1 - y 2 ) 2 . 

Here Xi and yi are the particle coordinates and a > 0 is the interaction 
strength. To define the self-adjoint Hamilton operator of the system 


we use the quadratic forms 


h- A m = 

/ VT 2 da;di/ and 


Jflxfl 

h v m = 

[ V\ty\ 2 dxd y, 


JClxQ. 

both defined on x H). 

Then by [2|, Theorem 1.8.1, the sum of 

the two forms has a closure h 

2 with 


h 2 f¥\ = /i_ A [*] + h v [ty] 

for all T in 

(1) Dom(/i 2 ) = Wq 2 {Q xfl)n Dom(hy). 

The positive self-adjoint operator associated with h 2 is 

H2 = -dl 1 -dl-dl 2 -d 2 y2 + V 

in L 2 (Vt x H). 

Lemma 2.2. 

a) For any choice of L , h and e one has 

a ess (H 2 ) C [V2a + 2n 2 , oo). 
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b) There is a choice of the constants L and h with L 2 + h 2 > 1 such 
that 

inf a(H 2 ) < \/2a + 27r 2 

for every e > 0, i.e., the operator H 2 has a bound state. 

From the above lemmata we conclude that a wave guide exists that 
has no bound state for one particle, but does have a geometrically 
induced bound state for two interacting particles. 

A remark on the physical interpretation of this effect is in order. 
As mentioned above, the argument of two particles acting like one of 
the double mass doesn’t apply to geometrically induced bound states, 
since their existence is mass-independent. To gain a physical intuition 
for our results anyway, we note that a bound state in a wave guide with 
bulges can be seen as a trade-off between reduced kinetic energy in the 
transverse direction (due to the enlarged cross-section) and increased 
kinetic energy in the longitudinal direction (due to the localization of 
the particle). Consider now two particles that attract each other and 
that would in free space form a ‘molecule’ with an average distance d 
between them. Assume for the case of our wave guide that d is con¬ 
siderably bigger than the cavity width h, but considerably smaller than 
the cavity length L. This means that in their transverse movement the 
two particles act rather as if they were independent of each other, thus 
receiving twice the energy decrease from the enlarged cross-section. In 
longitudinal direction, on the other hand, the two particles in the cav¬ 
ity behave like one particle of the double mass, such that the energy 
increase due to longitudinal localization is only half of what it would 
be for one particle alone. It follows that the energy trade-off is more 
‘favorable’ for the system of two interacting particles than for a single 
one. 


3. Two-particle bound states caused by Neumann 

BOUNDARY CONDITIONS 

If one introduces Neumann boundary conditions, an effect similar 
to the one described above happens even for particles in only one di¬ 
mension: Consider H 3 = —d 2 in L 2 (M + ) with a Neumann condition 
at x — 0. Then it is well known that cr ess (H 3 ) = M + and H 3 has no 
eigenvalues. Nevertheless, the corresponding two-particle Hamiltonian 
with an harmonic interaction turns out to have a bound state: 

We define the potential V = a\x\ — a^l 2 and the forms 
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on the restrictions of the functions in Qj°(M 2 ) to M + x M + . Then we 
can take h 4 to be the closure of h_ / \ + h V ] and its associated self-adjoint 
operator is 

h a = ~d 2 xi - d 2 X2 + a|xi - x 2 \ 2 

on L 2 (M + x M + ) with Neumann boundary conditions at x\ — 0 and at 
X 2 = 0 (see, e.g., [Sj, page 340). The domain of h 4 is 

(2) Dom (h 4 ) = W 1>2 (R + x M+) D Dom (h v ). 

Lemma 3.1. The operator H 4 has a bound state, i.e., an eigenvalue 
below the lower threshold of the essential spectrum. 

In view of Lemma Q it is no surprise that wave guides exist which 
have no one-particle bound states, but which do have a two-particle 
bound state induced by mixed boundary conditions. Omitting the 
proof, we only mention the simple example of a straight tube with 
Dirichlet boundary conditions on the edge and an additional Neumann 
condition imposed on one cross-section. 

4. Proofs of the results 

Proof of Lemma , \2.11 We introduce the operator Hi, which we define 
to be the Laplace operator on f2 with Dirichlet conditions on dfl and 
additional Neumann conditions on the set 

{(x, y) : x = ±L/2 and \y\ < |e|}. 

To prove Lemma 12.11 it is then sufficient to show that Hi has no spec¬ 
trum below 7T 2 . With the introduction of the new boundary conditions 
we have cut D into three separate domains: Two semi-strips and 
fl - in positive and negative x-direction, respectively, and the rectan¬ 
gle 0° = (—-§,-§) x (—t,|)- Thus H\ is the orthogonal sum of the 
Laplace operators on D + , fl - and D° (subject to appropriate boundary 
conditions), and 

a(Hi) = tr(-A n +) U tr(-A n -) U cr(-A n o). 

One can convince oneself easily that 

^(-An+) = o-(-An-) = [7 t 2 ,oo). 

The spectrum of — Aqo is purely discrete and if we call A(e) its lowest 
eigenvalue then 

A(0) = n 2 (h~ 2 + L~ 2 ) > 7T 2 . 

We can now apply a theorem of Gadyl’shin [Qj to see that A(e) — A(0) 
is of order e 2 , i.e., for small enough e > 0 we have inf cr(— Aqo) > n 2 . 
Altogether this means that 

inf cr(Hi) > inf a(H\) = if 2 for small e. 


□ 
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Proof of LemmaKE part a). Using the center of mass coordinates 

(3) u = ^ (x 2 + xi) and w = ^(x 2 - xi) 
we rewrite H 2 in the form 2 

(4) H 2 = ~^d 2 u - i<9 2 + 4 aw 2 - d 2 yx - d 2 y2 + a(y 2 - yf) 2 - 

To estimate the spectrum of H 2 from below we introduce Neumann 
boundary conditions on 

{(u,w,y!,y 2 ) : \w\ = (3} and {(u, w, y u y 2 ) : \w\ < (3, \u\ = (3 + ^}, 
for some (3 > 0, which turns H 2 into the orthogonal sum 

H 2 = H 2 \ {M>f3} © H 2 \ { ^ <f 3 ^ <fi+ L } © H 2 \ { ^ <( 3 M>/ 3 + l } . 

The spectrum of H 2 h, can be estimated from below by 4 a/3 2 

and the spectrum of H - 2 \,^ </3 ^ </ 3 + l^ is discrete. By separation of 
variables the spectrum of H 2 \n w ^ </3 ^ >f3+ Ly is found to be purely con¬ 
tinuous and its lower threshold is equal to the lowest eigenvalue of the 
‘transversal’ operator 

H t = ~dl + 4 aw 2 - d 2 yi - d 2 y2 + a(y 2 - yi) 2 

on L 2 ((—/3,/3) x (—1/2,1/2) 2 ) with Neumann conditions at \w\ = (3 
and Dirichlet conditions at \y\\ = 1/2 and \y 2 \ = 1/2. Neglecting the 
positive potential term a(y 2 — y±) 2 , we see that the lowest eigenvalue 
of H t is bigger than \p + 27T 2 , where \p is the lowest eigenvalue of the 
harmonic oscillator d ^ + 4 aw 2 on (—(3,/3) with Neumann boundary 
conditions. Below we will show that for (3 —> oo the eigenvalue A^ 
converges to a/2 a, i.e., the lowest eigenvalue of the harmonic oscillator 
on M. Consequently, for large enough (3 the lowest eigenvalue of H t is 
bigger than \[2a + 27T 2 . Part a) of Lemma 12.21 now follows from the 
fact that H 2 < II 2 and the min-max principle. 

It remains to show that linig^oo A^ = y/2a: Call hj = — i<9 2 + 4aw 2 
the Hamiltonian of the harmonic oscillator on the interval Id with 
Neumann boundary conditions. Then 

A /3 = inf a{h(-p t 0 )) = inf aih^ .-^©/^-/^©fy^oo)) < inf a(h R ) = V 2 a. 

The second step in the above chain of equalities follows from 

inf < 4a/3 2 and inf cr(h(_ 0O _ / g)) = inf cr(/i( j g j00 )) > 4 a/3 2 . 

Next we show that has a first eigenfunction that is sym¬ 

metric, non-negative and decreasing in |tc|: Let fp be a normalized 

2 In a slight abuse of notation we write H 2 for the two-particle Hamiltonian in 
Euclidean coordinates and for its unitarily equivalent counterpart in center of mass 
coordinates. 
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function such that h^p^^p = Xp4>p- We may assume that 4>p is ei¬ 
ther symmetric or antisymmetric, since otherwise we can replace it by 
4>p(w) + (f>p(— w). We write (j)*p for the symmetric decreasing rearrange¬ 
ment of <t>p (see m for the definition and properties of rearrange¬ 
ments). Then (f)*p is also normalized and belongs to the form domain 
W 1 ’ 2 ((-/3,/3))of Vftffl. The min-max principle yields 

r 0 

(5) 


Xp < 


< 


’^0 

r 

'-0 


+ ^W 2 (j)f^j dw 
QlVI 2 + 4 aw 2 0/^ dw = Xp. 


The second inequality in © follows from standard rearrangement the¬ 
orems 3 . The inequality is strict (and thus a contradiction) unless \(pp\ 
is decreasing in |w|. This shows that <j>p can be taken to be a non¬ 
negative symmetric eigenfunction to Xp that is decreasing in |w|. Then 
we have 


rp 


rP 


l-p 


4 aw 2 (j) 2 p{(5) dw < / 4 aw 2 (p 2 p(w) dw < Xp < vAha 


0 


and thus 
( 6 ) 

Now set 


0/3 M = 


MP) < 2- 5 / 4 3 1 / 2 a- 1 / 4 /?- 3 / 2 . 


<j> 0 (w) for \w\ < 0, 

(j)p(/3)(-\w\ + (3 + 1) for (3 < |w| < (3 + 1, 
0 for /3 + 1 < \w\ 


Then 


is in the form domain of /ir and we have 

Jr + daw 2 4 > 2 p(w] 


V 2 a = inf ct ( /i.r ) < 
rP + 1 

>0 


dw 


< A p + 2 


fu^Mdw 

\<i> 2 p{P) + 4 aw 2 (j) 2 p(f3) ) dw 


(7) — Xp + cj)p(f3) +-acj) 2 p((3)(3(3 2 + 3f3 + l) 

the penultimate step we used that f R (p 2 p(w) dw > f^p 4 >p{w) dw 


In 


= 1 


and the Ritz-Rayleigh characterization of A p. From © we conclude 


3 The estimate /(A, \4>p'\ 2 dw < /(A, \(j)p'\ 2 dw is a typical rearrangement property. 

It is usually stated for functions that go to zero at the boundary of their domain, but 
it also holds in the present case: Replacing </>p(w) by \<f>p(w)\ — \<f>p(j3)\ and 4>*p{w) by 

(\4>p (w) | — \<t>p(0)\)* does not change the value of the integrals, and \<f>p(w)\ — \<t>p(/3)\ 
is zero for w = ±/3 by (anti-) symmetry of 4>p. 
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that m converges to as (3 —> oo and therefore lim^oo \p = \[2a. 

□ 


Proof ofLemma\M V ar ^ b). We choose to fix the relations 
(8) a = L ~ 2 and hr 2 + L~ 2 M >1 


between the parameters that describe our wave guide. We define the 
domain hi as the set of all x 2 ,y 2 ) that satisfy the conditions 


u G 


(_3L 3 L\ 

V 8“’ “8“ ) ’ 



y 1 , 2/2 e 



J 


using the coordinates u and w as defined in (J3J). One can check that 
11 C SI x 11. We now define the test function T by 



on 0 and $ = 0 on (H x 0)\0, setting 

<t>(w) = e -^ w2 and C = <j>{L/ 8). 


Because the function T is Lipschitz continuous, has a bounded sup¬ 
port and vanishes at 5(0 x 11), we have G Wq 2 (Q x O). Since the 
potential V, restricted to the support of T, is bounded, we also have 
T G Dom(hy). By (JTJ) this means that T G Dom(/i 2 ). In the center of 
mass coordinates the quadratic form of H 2 reads 

h 2 m = I + (^ 1^) 2 + (< V ) 2 

+4crur|T| 2 + a(yi — 2 / 2 ) 2 1| J ) d w dudyi dy 2 . 


No we can apply the min-max principle with T as a test function to 
obtain 


inf a(H 2 ) < 



(9) 


+ 


87r 2 2 tt 2 

9L 2 + lP + 


7T 2 — 6 
67T 2 


ah 2 


I-L/siW^) 2 + 4 aw 2 ((j)(w) - C) 2 ) dw 
I-l/Mw) ~ C) 2 dw 


The last term can be estimated from above by 

I-l /+ 4 aw 2 cp(w) 2 ) dw 

!-L/M 2 i w ) ~ 2 C V>W) dw 
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which can, after an integration by parts, be written as 
V^a + (f-L/s dw ) 

l ~ 2 C (j-L%4?(w)dw) f^ 8 /s <f>(w)dw 

\J2 OL 

< -—- 

1 _ 26-^64 ^J_ L /8 g e -2V2 L-W jV8 g e ~V2L~^ dw 

where in the last step we have used that a = Za 2 and thus C = 
g-\/2i/64, Replacing w by the new variable w = w/y/L one can check 
that the product of the two integrals in the last line converges to a 
constant as L —> oo. Therefore, the last term in © can be estimated 
from above by a/2 ~a + 0(L~ 1 e~'^ L /® 4 ) for large enough L, which means 
that in view of (JHJ) 

inf a(H 2 ) <V2^+ 2Mvr 2 - ^ ^ + 0(L~ l e~^ Ll ^). 

If we choose L sufficiently large then the three last summands together 
are negative. If we then choose M sufficiently close to one, we get 
independently of e 

inf a(H 2 ) < V2a + 27r 2 , 

proving part b) of Lemma 12.21 □ 

Proof of Lemma \8.1\ In the center of mass coordinates H 4 acts in 
L 2 ({(u, w) : u > 0, |id| < u}) and takes the form 4 

H i =-\d 2 u~\dl + 4 aw 2 . 

Using a similar argument as in the proof of Lemma 12.21 part a), one 
can show that 

a ess (H A ) = [ V2a , oo). 

It remains to prove that H 4 has an eigenvalue below \/2a. We call 
4 >(w) the (positive and normalized) lowest eigenfunction of the har¬ 
monic oscillator — ^<9^+4 aw 2 in L 2 (M) and note that the corresponding 
eigenvalue is y/2a. We define the test function 

v L(u, w) = f>{w)e~ eu for u > 0, |u>| < u and some e > 0. 

We have T 6 W 1,2 ({(it, w) : u > 0, \w\ < u}) and since T drops off 
exponentially for u, |ru| —> oo, while V is only quadratic, also T e 


4 Again we abuse our notation and denote the two operators with respect to 
different coordinates by the same symbol H±, since they are unitarily equivalent. 







10 


HELMUT LINDE 


Dom(hy) holds. Thus T is in the form domain ((21) of H 4 and we can 
apply the min-max principle [T|j to obtain 

f u >o (|(<9 U T ) 2 + |(cCT ) 2 + 4 aw 2 #' 2 ) dtndw 

u > 0 T 2 dta d u 

|lt)| <14 

= 2 € + ' /2a + -/„>„ -' 

|iu| <1i 

I 11 the last step we used an integration by parts in w and the fact that 
0 satisfies the eigenvalue equation of the harmonic oscillator. The last 
summand is negative since <f{w) is positive, symmetric and decreasing 
in \w\, thus we have the estimate 


inf a(H 4 ) < 




/ 


infer ^4 < -e 2 + v / 2^+ W2 T \! * - 

v 2 f «>0 T 2 d w d u 

= -e 2 + + 2e f 0(u)0 / (-u)e _2e “ du 

2 J u > 0 

= V 2 a + e ( -e + 2 f 0 (u) 0 / (M)e _2e “ du J 
\2 A>o / 


The integral in the last line is negative and its absolute value increases 
when e goes to zero. Consequently, for some small enough e we have 
inf o(H/±) < y/2a, which proves Lemma fd.il □ 
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